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Summary. A singularly perturbed linear system of second order ordinary differen- 
tial equations of reaction-diffusion type with given boundary conditions is consid- 
ered. The leading term of each equation is multiplied by a small positive parameter. 
These parameters are assumed to be distinct. The components of the solution ex- 
hibit overlapping layers. Shishkin piecewise-uniform meshes are introduced, which 
are used in conjunction with a classical finite difference discretisation, to construct 
two numerical methods for solving this problem. It is proved that the numerical ap- 
proximations obtained with these methods are essentially first, respectively second, 
order convergent uniformly with respect to all of the parameters. 



1 Introduction 

The following two-point boundary value problem is considered for the singu- 
larly perturbed linear system of second order differential equations 

- Eu"{x) + A{x)u{x) ^i{x), x£{0,l), u(0) and u(l) given. (1) 

Here u is a column n — vector, E and A{x) are n x n matrices, 
E = diag(£), £ = (ei, • • • , £„) with < < 1 for all i = 1, . . . , n. The 
Si are assumed to be distinct and, for convenience, to have the ordering 

£1 < • • • < £„. 
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Cases with some of the parameters coincident are not considered here. 
The problem can also be written in the operator form 

Lu = f , u(0) and u(l) given 

where the operator L is defined by 

(f 

L = -ED^ + A(x) and D"^ ^ —. 

For all X G [0, 1] it is assumed that the components (x) of A(x) satisfy the 
inequalities 



aii{x) > ^|ay (a;)| for 1 < i < n, and aij(x) < for i 7^ j (2) 
and, for some a, 



j = l 



^ <<^ < (3) 

Wherever necessary the required smoothness of the problem data is assumed. 
It is also assumed, without loss of generality, that 

, \/a. , 

max 7^ < (4) 

The norms || V ||= maxi<fc<„ \Vk\ for any n-vector V, || y ||= supo<2,<i \y{x)\ 
for any scalar-valued function y and || y ||= maxi<fe<„ || yk || for any vector- 
valued function y are introduced. Throughout the paper C denotes a generic 
positive constant, which is independent of x and of all singular perturbation 
and discretization parameters. Furthermore, inequalities between vectors are 
understood in the componentwise sense. 



For a general introduction to parameter-uniform numerical methods for singu- 
lar perturbation problems, see [T], [5] and f?]. Parameter-uniform numerical 
methods for various special cases of ([Ij are examined in, for example, [5], 
[6] and [7]. For ([1]) itself parameter- uniform numerical methods of first and 
second order are considered in [8]. However, the present paper differs from 
[8] in two important ways. First of all, the meshes, and hence the numerical 
methods, used are different from those in [8]; the transition points between 
meshes of differing resolution are defined in a similar but different manner. 
The piecewise-uniform Shishkin meshes Mb in the present paper have the ele- 
gant property that they reduce to uniform meshes whenever b = 0. Secondly, 
the proofs of essentially first and second order parameter- uniform convergence 
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do not require the use of Green's function techniques, as is the case in [5]. The 
significance of this is that it is more hkely that such techniques can be ex- 
tended in future to problems in higher dimensions and to nonlinear problems, 
than is the case for proofs depending on Green's functions. It is also satisfy- 
ing, and appropriate in this special issue, to be able to demonstrate that the 
methods of proof pioneered by G. I. Shishkin can be extended successfully to 
problems of this kind. 

The plan of the paper is as follows. In the next section both standard and 
novel bounds on the smooth and singular components of the exact solution 
are obtained. The sharp estimates for the singular component in Lemma[5]are 
proved by mathematical induction, while an interesting ordering of the points 
established in Lemma |6l In Section 3 appropriate piecewise-uniform 
Shishkin meshes for essentially first order numerical methods are introduced, 
the discrete problem is defined and the discrete maximum principle and dis- 
crete stability properties are established. In Section 4 an expression for the 
local truncation error is found and two distinct standard estimates are stated. 
In Section 5 parameter-uniform estimates for the local truncation error of 
the smooth and singular components are obtained in a sequence of lemmas. 
The section culminates with the statement and proof of the essentially first 
order parameter-uniform error estimate. In the final section an outline of the 
construction and error estimation of an essentially second order parameter- 
uniform numerical method is presented. 



2 Analytical results 

The operator L satisfies the following maximum principle 

Lemma 1. Let A{x) satisfy (0j and Let ■0 be any function in the 
domain of L such that ^p{0) > and ■0(1) > 0. Then Lxp{x) > for all 
X e (0,1) implies that xp{x) >0 for all x E [0,1]. 

Proof. Let i*,x* be such that 'ipi*{x*) = min^^a, -0^ (x) and assume that the 
lemma is false. Then 0^. (a;*) < . From the hypotheses we have x* ^ {0, 1} 
and 0-'.(a;*) > 0. Thus 

n 

(L0(a;*))i. = -ei*i/l,{x*) + ^a,.,j{x*)ipj{x*) < 0, 

which contradicts the assumption and proves the result for L. ■ 

Let A{x) be any principal sub-matrix of A{x) and L the corresponding op- 
erator. To see that any L satisfies the same maximum principle as L, it sufiices 
to observe that the elements of A{x) satisfy a fortiori the same inequalities 
as those of A{x). 
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We remark that the maximum principle is not necessary for the results 
that follow, but it is a convenient tool in their proof. 

Lemma 2. Let A{x) satisfy (0) and If "^p is any function in the domain 
of L, then for each i, I < i < n, 

\^P,{x)\ < maxjll ^(0) ||,|| |1 L^P [0,1]- 

Proof. Define the two functions 

= maxjll ^P{0) ||, || || , i || || | e ± ,p{x) 

where e — (1, 1)-'" is the unit column vector. Using the properties of 
A it is not hard to verify that 9^(0) > 0, > and L0±(x) > 0. 

It follows from Lemma [T] that 6^{x) > for all x e [0, 1]. ■ 

A standard estimate of the exact solution and its derivatives is contained in 
the following lemma. 

Lemma 3. Let A(x) satisfy (0j and and let u be the exact solution of 
Then, for each i = I ... n, all x Cz [0,1] and fc = 0, 1, 2, 

|"f (x)|<Cert(||u(0)|| + ||u(l)|| + ||f||) 

and 3 

\u'"(x)\ < Ce~n||u(0)|| + ||u(l)|| + llfll + llf'll) 



Proof. The bound on u is an immediate consequence of Lemma [2] and the 
differential equation. 

To bound u'^{x), for all i and any x, consider an interval Nx — [a,a + y/ei] 
such that X G Nx. Then, by the mean value theorem, for some y g Nx, 

I Ui{a + Jel) - Ui{a) 
u^{y) = ^-—= 

and it follows that 

K(2/)|<2£-'|k||. 

Now 

u'{x)^u'{y)+ u"{s)ds = u'{y) + E-^ {-i{s) + A{s)u{s))ds 
Jy Jy 

and so 

\<{x)\ < \u'M\ + Ce-\m\ + Hull) fds < Ce-H\\f,\\ + ||u||) 



from which the required bound follows. 

Rewriting and differentiating the differential equation gives u" — E~^{Au~i), 
u'" = E-^{Au' + A'u - {'), and the bounds on u'/, uf follow. ■ 
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The reduced solution Uq of H]) is the solution of the reduced equation 
Auq = f. The Shishkin decomposition of the exact solution u of ([1} is 
u = V + w where the smooth component v is the solution of Lv — f 
with v(0) — Uo(0) and v(l) — Uo(l) and the singular component w is 
the solution of Lw = with w(0) = u(0) - v(0) and w(l) = u(l) - v(l). 
For convenience the left and right boundary layers of w are separated using 
the further decomposition w = w' + w'' where Lw' ~ 0, w'(0) = u(0) — 
v(0), w'(l) = and Lw'' = 0, w''(0) = 0, w''(l) = u(l) - v(l). 
Bounds on the smooth component and its derivatives are contained in 

Lemma 4. Let A{x) satisfy 0^ and Then the smooth component v and 
its derivatives satisfy, for all x ^ [0,1], i — 1, . . . n and k = 0, ... 3, 

\vl'\x)\<Cil+s]-^). 

Proof. The bound on v is an immediate consequence of the defining equations 
for V and Lemma [H 

The bounds on v' and v" arc found as follows. Differentiating twice the equa- 
tion for V, it is not hard to see that v" satisfies 

Lv" = g, where g = f " - A"v - 2A'v'- 

Also the defining equations for v yield at a; = 0, x — 1 

v"(0) = O, v"(l) = 0. 

Applying Lemma [2] to v" then gives 

||v"||<C(l + ||v'||). (5) 

Choosing i*, x* , such that I < i* < n, x* ^ (0, 1) and 

^^(^*) = l|v'|| (6) 

and using a Taylor expansion it follows that, for any y G [0, 1 — a;*] and some 
r], X* < rj < X* + y, 

Vi-{x* +v) = v,.{x*) + y^^'{x*) + y ('7)- (7) 
Rearranging ([7]) yields 

/ , *^ V^' {X* + y) ~ V,> (X*) y „ , . 

v^-{x ) = 7;v,,{v) (8) 

y ^ 



and so, from ([6]) and 



<-||v|| + ^||v"||. (9) 

y 2 
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Using dni), ^ and the bound on v yields 



(l-^)||v"||<C(l + -). 

2 y 



(10) 



Choosing y = min(^, 1-a;*), (UHl) then gives | |v"| | < C and (HD gives ||v'|| < C 
as required. The bound on v'" is obtained by a similar argument. ■ 



The layer functions S', , i — 1, 
are defined on [0, 1] by 



associated with the solution u, 



Blix)^e-^V^^, Bl{x) = B\{l-x). 

The following elementary properties of these layer functions, for all 1 < j < 
j < n and < x < y < 1, should be noted: 

(i) Blix) < B]{x), B\{x) > Blivl < Blix) < 1. 

(ii) Blix) < Blix) < Bliy), < < 1. 

Bounds on the singular components w' , w'' of u and their derivatives are 
contained in 

Lemma 5. Let Aix) satisfy (0) and ^.Then there exists a constant C, such 
that, for each x G [0, 1] and i — 1, . . . , n, 



\wlix)\ < CBUx), w\''ix) 



< C 



B[ix) 



q=i 



Wi ix) 



<c±m. 



SiW. 



1,111 



ix) 



q=l 



Analogous results hold for wl and its derivatives. 

Proof. First we obtain the bound on w'. We define the two functions 6^ = 
CBle±yvK Then clearly 0±(O) > 0, 6/±(l) > and L6»± ^ CLiB^e). Then, 

)Bl > 0. By Lemma [II 6^ > 0, 



for i = 1, . . . , n, iLe^)i = C(X;"=i a^j - aj 
which leads to the required bound on w' 

Assuming, for the moment, the bounds on the first derivatives w'"', the 
system of differential equations satisfied by w' is differentiated to get 



-Ew 



Llll 



Aw 



I,' 



A'w' = 0. 



The required bounds on the w''"' follow from those on and w^'' . It remains 
therefore to establish the bounds on w\'' and w\'" , for which the following 
mathematical induction argument is used. It is assumed that the bounds hold 
for all systems up to order n — 1. It is then shown that the bounds hold for 
order n. The induction argument is completed by observing that the bounds 
for the scalar case n = 1 are proved in [1]. 



I,' 
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It is now shown that under the induction hypothesis the required bounds 
hold for w\' and w\" . The bounds when i=n are established first. The differ- 
ential equation for w\ gives Snwli" = (Aw')„ and the required bound on wli" 
follows at once from that for w'. For wli' it is seen from the bounds in Lemma 

m applied to the system satisfied by w', that [^'''(a;)! < Ce.^ ^ . In particu- 

lar, |w5j'(0)| < Cen ^ and < Cen ^ . It is also not hard to verify that 

Lw''' = — A'w'. Using these results, the inequalities i < n, and the 

properties of A, it follows that the two barrier functions 6^ — Ci?^ 2 i?^e±w''' 
satisfy the inequalities 9^(0) > 0, > 0, and L0± > 0. It follows from 

Lemma [T] that 6^ > and in particular that its n*'' component satisfies 

_i 

< e'en ^-B^(x) as required. 
To bound wl'' and for 1 < i < n — 1 introduce w' = {w{, . . . ,'wl^_i). 
Then, taking the first n ~ 1 equations satisfied by w', it follows that 

-E'w'^" + iw' = g, 

where E, A is the matrix obtained by deleting the last row and column 
from E, A, respectively, and the components of g are gi = —ai^nw\^ for 
1 < i < n — 1. Using the bounds already obtained for wl^,wli', it is seen 

that g is bounded by CB\^{x) and g' by C^^^. The boundary conditions 

for w' are w' (0) = u(0) - u°(0), w'(l) = 0, where is the solution of the 
reduced problem u° = A^^i, and are bounded by C(|| u(0) || + || f(0) ||) 
and C(|| u(l) || + || f(l) ||). Now decompose w' into smooth and singular 
components to get 

w' = q + r, w'-' = q' + r'. 

Applying Lemma 3] to q and using the bounds on the inhomogeneous term g 
and its derivative g', it follows that |q'(a;)| < and |q"(a;)| < C^^. 

Using mathematical induction, assume that the result holds for all systems 
with n — 1 equations. Then Lemma [5] applies to r and so, for i ~ l,...,7i — 1, 

K{x)\<cY.^, K'(.)|<cg4^. 

Combining the bounds for the derivatives of qi and r^, it fohows that 

q—l V ^ q—i ^ 

Thus, the bounds on w\'' and w''/' hold for a system with n equations, as 

required. The proof of the analogous results for the right boundary layer 
functions is analogous. ■ 

Definition 1. For B\, Bj and each 1 < i ^ j < n, the point defined by 
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It is remarked that 



(11) 



In the next lemma it is shown that the points Xi^j exist, are uniquely 
defined, lie in the domain [0, i] and have an interesting ordering. 

Lemma 6. Assume that holds. If, in addition, y/el < then, for all 

i,j with 1 < i < j < n, the points Xij exist, are uniquely defined, lie in (0, ^] 
and satisfy the following inequalities 



In addition the following ordering holds 

Xij < Xi^ij, if i + 1 < j and Xij < x^.j+i, if i < j. (14) 



Analogous results hold for the , Bj and the points 1 — Xij. 

Proof. Existence, uniqueness and (|13p follow from the observation that < 
y/Sj, for i < j, and the ratio of the two sides of (fTTj) . namely 

' — T— — = exp { — \/ax{^— ), 

is monotonically decreasing from the value -^S > 1 as x increases from 0. 
The point the unique point x at which this ratio has the value 1. 

Rearranging pTj) gives 



^) ^) 



(15) 



Using the hypotheses it follows that 



x.,<^ln(^)<Ml^^.!^^<i 

as required. 

To prove returning to and writing = exp(— p^), for some pk > 
and all fc, it follows that 
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The inequality Xij < Xi+ij is equivalent to 

Pi - P] ^ Pt+i - Pj 
exppi-exppj exppi+i - exppj ' 

which can be written in the form 

{Pi+l - P]) exp(p, - Pj) + {p^ - Pi+i) - {p^ - Pj)exp{p^+i - Pj) > 0. 

With a — Pi — Pj and b = Pi+i — Pj it is not hard to see that a > 6 > and 
a — b — Pi — Pi+i- Moreover, the previous inequality is then equivalent to 

exp a — 1 exp b — \ 
a b ' 

which is true because a> b and proves the first part of (I14p . The second part 
is proved by a similar argument. 

The analogous results for the -B[, and the points 1 — cc^j are proved by a 
similar argument. 



3 The discrete problem 

A piecewise uniform mesh with N mesh-intervals and mesh-points {I'ij^o 
is now constructed by dividing the interval [0, 1] into 2n + I sub-intervals as 
follows 

[0, di] U • • • U (cr„_i, a„] U (a„, 1 - cr„] U (1 ~ (7„, 1 - cr„_i] U • • • U (1 - <Ti, 1]. 

The n transition parameters, which determine the points separating the uni- 
form meshes, are defined by 



cr„ = mm i, 



(16) 



and for i — 1, 



, n — 1 



cr, = mm 



O'i+l 



a 



(17) 



Clearly 



< (Tl < 



i 

< o-„ < -, 



^ < 1 - cr„ < 



1 - (Tl < 1. 



Then, on the sub- interval ((t„,1 — cr„] a uniform mesh with mesh- intervals 
is placed, on each of the sub-intervals (ctj, ctj+i] and (1 — Ci+i, 1 — ai], i — 
I, ... ,n — 1, a uniform mesh of 2"-^+'^ mesh-intervals is placed and on 
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both of the sub-intervals [0, CTi] and (1 — cti, 1] a uniform mesh of -^^^ 
mesh- intervals is placed. In practice it is convenient to take N = 2"-k where k 
is some positive power of 2. This construction leads to a class of 2" piecewise 
uniform Shishkin meshes Mb, where b denotes an n-vector with 6^ = if 
Ci = and hi = 1 otherwise. Note that Mb is a classical uniform mesh 
when b = 0. It is not hard to see also that 

a^ = 2"(^~'+iVj+i when h, = --- = hj=0. (18) 

B\{(Ji) = Bl{l - ai) = N-^ when b^ = 1. (19) 
Writing 5j = Xj+i — Xj-i note that, on any Mb, 

5j<CN~^, l<j<N^\ (20) 

and 

o-i < C^ellniV, l<i<n. (21) 

Furthermore, 

B\{ai — {xj — Xj-i)) < CB\{ai) if Xj — cji, for some (22) 

To verify note that if xj — <Ji then Xj — Xj^i = (j^/2''-i~+^) — "'^o'iS"^*^^ 
and the result follows from Bl{ai — (xj — Xj-i)) < B\{ai — < 
exp(i|lli^)i?Ka,) and exp(i|^iii^) < C. 

The discrete two-point boundary value problem is now defined on any mesh 
Mb by the finite difference method 

~ES^lJ + A{x)V:^i{x), U(0) =u(0), U(l) = u(l). (23) 

This is used to compute numerical approximations to the exact solution of 
([T]). Note that ([23|l can also be written in the operator form 

L^U = f, U(0) = u(0), U(l)=u(l) 

where 

= ^Ed^ + A{x) 

and S^, £)+ and arc the difference operators 



S-'Vix 



D+U{xj)-D-\J{x,) 



^+U(.,) = and D-lJix,) = U(^.)-U(^.-i). 

hj+i hj 

with h-j = , h.j = Xj-Xj-i. 

The following discrete results are analogous to those for the continuous case. 
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Lemma 7. Let A{x) satisfy (0j and Q). Then, for any mesh function the 
inequalities "^{0) > 0, *(1) > OandL'^*(xj) > /or 1 < j < - 1, 
imply that ^(a;^) > /or < j < iV. 

Proof. Let i*,j* be such that ^i*(xj*) = mhii j^i{xj) and assume that the 
lemma is false. Then ^i*{xj*) < . From the hypotheses we have j* ^ 0, N 
eind'l',.{xj.)-<I'^.{xj.-i) <0, > 0, so 6^W„{xj^) > 0. 

It follows that 

n 

(L''^*(a;j*)) = -ei*(5^tf'i*(2;j*) + ^Oi,, fc(a;j,)!f'fc(a;j*) < 0, 

k=l 

which is a contradiction, as required. 

An immediate consequence of this is the following discrete stability result. 

Lemma 8. Let A{x) satisfy and (0). Then, for any mesh function , 

II ^{x,) II < max|||*(0)|U|*(l)||, ^l|L^*ll} , < j < N. 

Proof. Define the two functions 

0±(x,) = max{||*(0)||, ||*(1)||, i||LN*||}e± 

a 

where e = (1, ... , 1) is the unit vector. Using the properties of A it is not 
hard to verify that 0±(O) > 0, 0=^(1) > and L^@^{x.j) > 0. It follows 
from Lemma [7] that 0=^(xj ) > for all < j < iV. 

4 The local truncation error 

From LemmalHl it is seen that in order to bound the error | |U — u| | it suffices 
to bound L^(U — u). But this expression satisfies 

L^(U - u) = L^(U) - L^(u) = f - L^(u) = L(u) - L^(u) 

= (L - L^)u = -E{S^ - D^)u 
which is the local truncation of the second derivative. Then 

E{d^-D^)u = E{S^ - D^)v + E{S^ - D^)w 

and so, by the triangle inequality, 

II L^(U - u) II < II E{S' - D')^ II + II E{6' - D')^ \\ . (24) 
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Thus, the smooth and singular components of the local truncation error can be 
treated separately. In view of this it is noted that, for any smooth function ip, 
the following two distinct estimates of the local truncation error of its second 
derivative hold 

1(^2 _ D^)^{xj)\ < 2 max \^"{s)\ (25) 
s e Ij 

and 

US' - D^)ij{x,)\ < ^niax|V'"'(s)| (26) 
o s e Ij 

where Ij — [xj^i,Xj+i]. 



5 Error estimate 

The smooth component of the local truncation error is estimated in the fol- 
lowing lemma. 

Lemma 9. Let A(x) satisfy 0) and (0). Then, for each i — I, . . . , n and 
j — 1, . . . , N — 1, we have 

\e,iS' ~ D')v,{xj)\ < C^^N-\ 

Proof. Using ([261), Lemma Hand it follows that 

\e^{5'' - D'')v,{x,)\ < C<5j max |£,<'(s)| < C^,5j < C^^N-^ 

s e Ij 

as required. ■ 

For the singular component a similar estimate is needed, but in the proof 
the different types of mesh must be distinguished. The following preliminary 
lemmas are required. 

Lemma 10. Let A{x) satisfy 0) and (0). Then, for each i = 1, . . . , n and 

j — 1, . . . , N , on each mesh M^, the following estimate holds 

\e.iS' - D')wlix,)\ < ^. 

Proof. From (|26p and Lemma [5l it follows that 

\e,iS' - D')wl{xj)\ < C6j max \e,w/" {s)\ 

S G Ij 



< CS j max > V 

C6, 



< 

as required. ■ 
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In what follows second degree polynomials of the form 

ViAx) = 2^ [xe) 

k=0 

are used, where 6 denotes a pair of integers separated by a comma. 

Lemma 11. Let A{x) satisfy ^ and Then, for each i = 1, n, 
J = 1, . . . , TV and k = 1, . . . , n — 1, on each mesh Mb with bk — 1, there 
exists a decomposition 

k+l 
m— 1 

for which the following estimates hold for each m, 1 < m < k, 
\e.w'Ux)\ < CBUx), \e.w'lUx)\ < C^''^ 

V ^' 

and 

<;=fe+l 



Furthermore 

Analogous results hold for the and their derivatives. 

Proof. Since 6^ = 1 it follows that y/Fk < y/ek+i/2, so Xk,k+i G (0, ^) and 
the decomposition 

fc+i 

m— 1 

exists, where the components of the decomposition are defined by 

Pi-k,k+i on [0,Xk,k+i) 



^^■'=+1 ~ 1 otherwise 



and for each m, k > m > 2, 

on [0,X,„_i,,n) 

- Egii+i Otherwise 
and 

fc+i 

Wi,i ^ wl -'^Wi^q on [0,1]. 

g=2 
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From the above definitions it follows that, for each m, I < m < k, Wi 

on 

To establish the bounds on the third derivatives it is seen that: 
for X G [xk,k+i, 1], Lemma [S] and x > Xk,k+i imply that 



k.<Ui(-)l = \e.wr'{x)\ < ^ < C ^ 



9=1 q=k+l ^ 



for X e [0,a;fe_fc+i], Lemma [S] and x < Xk^k+i imply that 

= k^U;,' (a^fc,fc+l)| < 2^ ' < 2^ _ < 2^ 

and for each m = k, ... , 2, it follows that 

for X e [a;m,m+i, 1], = 0; 

for X e [xrn-i.m, Xm,m+i], Lemma [5] implies that 



for X e [0,a;m-i,m]i Lemma [5] and x < Xm_i .jn imply that 

I /// / M I 1,11'/ ^, ^ ^Sr^ Bli^rn~l,m) B\ {Xjn-l,m) ^^BI^{x) 
= \£tw- < — < C — < C 



for X e [a;i,2, 1], w-" = 0; 

for X S [0,a;i^2], Lemma [5] implies that 



fc+l n 



k.<'i(x)| < \e.wY"{x)\ + ^ |e.<;(x)| < ^ < C 
For the bounds on the second derivatives note that, for each m, 1 < m < k 

for X e [Xrn,m+1, 1], = 0; 

for X e [0,a:TO,,„+i], J^^^-^+i e,w^"^(s)ds = eiW^^„(x™^„i+i)-£iW-'„(a;) = 
-Siw'l^ix) 
and so 

< r"^'"^^' \e.w'/'{s)\ds < BUs)ds < CBI{x). 



Finally, since 

k 

\e,{5^ - D^)w\{xj)\ < J2 - D^)w^,m{xJ)\ + \e,{S^ - D^)w,,k+i{xj)\, 
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using (pS)) on the last term and on all other terms on the right hand side, 
it follows that 

k 

|e,(52 _ D'^)w\{xj)\ < CiY, max|e,u;^;„^(s)| +^,max|e,<'fe+i(s)|). 

m—l 

The desired result follows by applying the bounds on the derivatives in the 
first part of this lemma. 

The proof for the and their derivatives is similar. 

Lemma 12. Let A(x) satisfy 0) and (0). Then, for each i — 1, . . . , n and 
j = 1, . . . , N , on each mesh Mh the following estimate holds 

\e,{S^-D^)wl{x,)\<CBl{x,^^). 

An analogous result holds for . 

Proof. From ((25l) and Lemma [51 for each i = 1 , . . . , n and j = 1, . . . ,N, 
it follows that 

\e,{S^ - D^)wl{xj)\ < C max|£,«;^^"(s)| 

< Ce^Y.^^^^^ ^ CBUx,^,).. 

p—t ^ 

Using the above preliminary lemmas on appropriate subintervals, the de- 
sired estimate of the singular components of the local truncation error are 
proved in the following lemma. 

Lemma 13. Let A{x) satisfy 0j and (0). Then, for each i ~ 1, . . . , n and 
J = 1, . . . , N , the following estimate holds 

Isiid^ - D^)wi{xj)\ < CN^^ InN. 

Proof. Since w = w' + w'', it suffices to prove the result for w' and w'' 
separately. Here it is proved for w' . A similar proof holds for w'' . 
Stepping out from the origin each subinterval is treated separately. 
First, consider x £ (0, (Ti). Then, on each mesh Afb, Sj < CN~-^ai and the 
result follows from (|2ip and Lemma [TUl 

Secondly, consider x e {(Ji,(J2), then (Ti < xj^i and 6j < CN^^a2. The 2"+^ 
possible meshes are divided into 2 subclasses. On the meshes Mb with 6i = 
the result follows from ([2T|l . (fT8| and Lemma [TOl On the meshes Mb with 
bi — 1 the result follows from ([2T|) . (fT9|) and Lemma [TTl When x — ai, similar 
arguments apply for the 2 subclasses, except that ([22|) is also needed for the 
second subclass. 

Thirdly, in the general case x € (<7m,(7m+i) for 2 < m < n — I, it follows 
that (T,„ < Xj_i and Sj < CN~^a,n+i- Then Mb is divided into 3 subclasses: 
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- {Afb : foi = • • • = b,„ = 0}, Af£ = {A/b : K ^ 1, K+i - • • • = 
6,„ = for some 1 < r < to - 1} and M^' = {Mb : &m = !}• On A/° the 
result follows from (HIl), (HE]) and Lemma [TUl on from (PT|) . P^ . (fTO)) and 
Lemma [11] on from (pi]) , (fT5|) and Lemma [TlJ When a; = (7™ , similar 
arguments apply for the 3 subclasses, except that (l22l) is also needed for the 
third subclass. 

Finally, for x £ {an, 1), cr„ < Xj-i and 6j < CN^^. Then Afb is divided into 
3 subclasses: A^^ = {Afb : &i = • • • = = 0}, A/^; = {Afb : &r = 1, ^'r+i = 
••• = &„ = for some 1 < r < n - 1} and = {Afb : 6„ = 1}. On 
the result follows from ([21]), (HH]) and Lemma [H on Af£ from ((2T1) . pH) . (fTg]) 
and Lemma [TTl on Afb from p9|) and Lemma fT2l When x = cr„, similar 
arguments apply for the 3 subclasses, except that ((22|) is also needed for the 
third subclass. ■ 

Let u denote the exact solution from ([T]) and U the discrete solution from 
(|23p . Then, the main result of this paper is the following parameter uniform 
error estimate 

Theorem 1. Let A{x) satisfy 0) and Q). Then there exists a constant C 
such that 

II U - u ||< CiV^MniV, 

for all N >l. 

Proof. This follows immediately by applying Lemmas [5] and [T3] to ([M]) and 
using Lemma [S] ■ 



6 An essentially second order method 

In this section it is shown that a simple modification to the Shishkin mesh 
constructed above leads to an essentially second order parameter-uniform nu- 
merical method for ^ . The finite difference operator is the same as for the first 
order method; the Shishkin piecewise uniform mesh is modified by choosing 
a different set of transition parameters. Instead of (fT6|) and (fT7|) the following 
parameters are used 




and for i = 1 , ... , n — 1 

T,=min|^,2y^lniv|. (28) 

The proof that the resulting numerical method is essentially second order 
parameter-uniform is similar to the above and is based on an extension of the 
techniques employed in [3]. It is assumed henceforth that the problem data 
satisfy additional smoothness conditions, as required. 
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It is first noted that, with the definitions (P7| . (pS)) . is replaced by 

Blin) = Hr(l - n) = when 6, - 1. (29) 

Also, for any smooth function tp, in addition to ([^5]) and (pS)) . the following 
estimate of the local truncation error holds at any mesh point Xj of a locally 
uniform mesh 

\{S^ - D^)iP{xj)\ < max \ip""is)\ if a;^ - Xj-i = xj+i - Xj. (30) 
12 s e /j 

From their construction, it is clear that the above Shishkin meshes are locally 
uniform everywhere, except at the points Xj = where k € I^, and /b = {fc : 
bk = 1}. 

To estimate the smooth component of the error, note that the estimate in 
Lemma [9] can be modified to 

Now introduce the mesh function # where, for each I < i < n, 

<P,{x,) = CN-^{e{xj) + 1), 
where 9 = X]fee/b ^'^ a-nd, for fc e /b, 9k is the piecewise constant polynomial 

r 0, X e [0, Tfe) 

6'fe(x) = < 1, X e [Tfe, 1 - Tfc] 

[ 0, (1-Tfc,i] 

Then 

< <P^{xj) < CN~^, l<i<n 

and 

n 

(L^^>(x,)). = CN-'[-e,6^9{x,) + Y,a^A^,m^,) + !)]■ 
It follows that 

Considering the cases Si > N~-^ and Si < separately, choosing C" suffi- 
ciently small, comparing (j32p with (|3ip and applying the discrete maximum 
principle to the barrier functions 

*± = * ± (V - v) 

gives the following estimate 
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\\(V-^)\\<CN-^. (33) 

To estimate the singular component of the error, note that the estimates 
in Lemmas [10] and [11] are modified, respectively, to 

\e,{S' - D')wi{x,)\ < (34) 
\e.{S^ - D^)wlix,)\ < C{Bi{x,^,) + (35) 

£k+l 

Combining these with Lemma [T2l and repeating the same for the wl leads to 
the following estimate of the singular component of the local truncation error 

|e,(^2 _ D'^)w^{xj)\ < C{N-^ \nN)^. (36) 

Application of Lemma [8] then gives 

||(W - w)|| < C(iV~MniV)^ (37) 

Combining (|33p and (|37|) leads at once to the required essentially second order 
parameter-uniform error estimate 

\\{V -u)\\<CiN-HnNf. (38) 
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